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Adsorption of a symmetric (AB) random copolymer (RC) onto a symmetric (ab) random heterogeneous surface 
(RS) is studied in the annealed approximation by using a two-dimensional partially directed walk model of the 
polymer. We show that in the symmetric case, the expected a posteriori compositions of the RC and the RS 
have correct values (corresponding to their a priori probabilities) and do not change with the temperature, 
whereas second moments of monomers and sites distributions in the RC and RS change. This indicates that 
monomers and sites do not interconvert but only rearrange in order to provide better matching between them 
and, as a result, a stronger adsorption of the RC on the RS. However, any violation of the system symmetry 
shifts equilibrium towards the major component and/or more favorable contacts and leads to interconversion 
of monomers and sites. 
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1. Introduction 

Adsorption of heteropolymers — polymers composed of monomer units of two and more types — 
onto chemically heterogeneous surfaces was intensively studied in the last decades. A particular interest 
to this problem is motivated by its connection to the question of molecular recognition playing a cru¬ 
cial role in living organisms and in various biomedical/biotechnological applications. To understand the 
mechanisms of the polymer-surface recognition, the problem was extensively investigated from different 
angles by using relatively simple and physically transparent models. 

In particular, directed walk models of polymers (l}{^ played an important role in studying homopoly¬ 
mer adsorption 14HT21 and coUapse 11311171 and related problems of force-induced desorption I18H26I and 
unfolding of a collapsed macromolecule (27l|^. Random 12911311 and periodic 13211341 copolymer adsorp¬ 
tion and mechanical desorption 13511371 were also studied with the aid of directed models. The advantage 
of directed models consists in their simplicity; polymer directedness allows one to obtain an exact solu¬ 
tion in most cases, typically in the long chain limit. At the same time, they provide a physically reasonable 
and tractable picture of the phenomenon under study which is commonly in agreement with the results 
of a more realistic computer simulation of the same (or similar) system in 3 dimensions. Another attrac¬ 
tive feature of directed polymer models is the inherent self-avoidance of polymer conformations (and the 
corresponding impossibility of visiting the same surface site by two different monomer units simultane¬ 
ously). 

To study homo- or heteropolymer adsorption onto a homogeneous surface, one can use fully directed 
polymer models (in particular, Dyck or Motzkin paths). In the case when both the polymer and the sur¬ 
face are heterogeneous and the polymer should adjust its conformation to the surface pattern in order 
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to attain a better matching of its monomer sequence to the heterogeneous surface pattern and maximize 
the amount of favorable polymer-surface contacts, the minimal model that allows one to consider this 
phenomenon is the two-dimensional partially directed walk (2D-PDW) model. The fully directed model, 
which is simpler than the 2D-PDW one, is not suitable for this purpose because for a particular monomer 
unit, there is one and only one surface site that it can visit (hence, heteropolymer adsorption onto het¬ 
erogeneous surface is equivalent to the situation of a heteropolymer adsorption onto a homogeneous 
surface). 

In our recent paper (HD, a random copolymer (RC) adsorption onto a random surface (RS) was consid¬ 
ered in the framework of the 2D-PDW model of the polymer on a square lattice (the “surface”, therefore, 
was simply a line). In order to take correlations into account, both random sequences of monomers (in 
the RC) and sites (in the RS) were modelled as first-order Markov chains. The problem was solved by using 
a combination of the annealed approximation to perform double averaging over sequence and surface 
disorder, and the generating functions (GFs) approach to sum over aU conformations of the RC. The key 
result of the work was the derivation of an equation to find the smallest singularity of the GF of 
the adsorbed chain. The latter provides an asymptotic form of the canonical partition function for the 
annealed system which, in turn, gives an access to the calculation of various observables. This equation 
was then apphed in for the analysis of the adsorption transition point for different sets of the system 

parameters. This allowed us to study the effect of the interplay between correlations in the RC and the RS 
on the transition temperature. 

In the present work, we employ the model introduced in (31) for a comprehensive study of the RC 
adsorption onto the RS beyond the transition point in the annealed approximation. It is important to note 
that for the system considered, the annealed approximation is interesting not only as a mathematical trick 
but due to its correspondence to real physical situations. Mathematically, the annealed approximation 
is equivalent to direct averaging the partition function of a disordered system instead of averaging its 
logarithm (i.e., free energy) over all possible realizations of the disorder. For our system this means that 
monomers and sites participate in thermal motion along with the conformational degrees of freedom, 
and, following Grosberg (38), we may refer to the annealed RC (RS) as a copolymer (a heterogeneous 
surface) with a “mobile primary structure”. 

Yoshinaga et al. (33 developed a theory of adsorption of the so-called two-state polymers consisting of 
monomer units that can change its type (i.e., its affinity to the substrate or its hydrophobicity/hydrophili- 
city). The authors of (39) showed the equivalence of the two-state polymer to the annealed two-letter RC 
and established the correspondence between the RC parameters (a priori probabilities to find a monomer 
of a certain type in the RC sequence) and standard chemical potentials of monomers in the bulk and at the 
surface. Another situation that can be described in terms of the annealed approximation was suggested 
in (38): the annealed RC can be viewed as a homopolymer consisting of monomer units that can adsorb 
smaU molecules, like surfactants gO), from the solution. As a result, there are two types of monomer 
units: free and with an adsorbed “side group”. Correspondingly, monomer units with “side groups” may 
have additional attraction to each other or to a surface. This picture can be straightforwardly extended to 
the case where two interacting objects are heterogeneous: both polymer and surface may be “two state” 
or may bind (different) ligands and, being in a bound state, attract each other. 

There exists another aspect that makes the study of RC adsorption onto RS in the annealed approxi¬ 
mation interesting. As it is weU known, in the case of a RC adsorbing onto a homogeneous surface, the 
situation develops as foUows (29)141): with a decreasing temperature, transformation of non-adsorbing 
monomer units into adsorbing ones occurs. As a result, the a posteriori first and higher moments of the 
monomer distribution (that is, the number of monomer units of A and B type, the number of dyads: AA, 
BB, AB and BA, triads: AAA, ABA, BAA,... and so on) do not correspond to their expectation values. In the 
case where both interacting objects are heterogeneous and there are two or more possibilities of forming 
favorable, or “good”, contacts (say, there are two kinds of favorable contacts: Aa and 5b) it is not easy to 
predict the system behavior. The most unclear will be the symmetric situation, where both good contacts 
(Aa and Bb) are equally favorable, whereas non-attractive “bad” contacts (Ab and Ba) are equally unfavor¬ 
able and, in addition, both the RC and the RS have a symmetric composition (i.e., equal amounts of A and 
B monomers and a and b sites): here, one cannot say in advance how this compositional equilibrium wiU 
be biased during the interaction of the RC with the RS upon a decreasing temperature (increasing inter¬ 
action strength). Phase diagram for the symmetric case was analyzed in (31): it was shown that RC tends 
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to adsorb oitto the RS with the same type of correlations: quasi-blocky on quasi-blocky, quasi-alternating 
on quasi-alternating whereas uncorrelated, or Bernoullian, RCs (RSs) do not “feel” correlations on the RS 
(in the RC). These findings are in line with the results of Monte Carlo simulations (42l for a more realistic 
three-dimensional self-avoiding chain. 

The rest of the paper is organized as follows. In section]^ the model and the formalism based on the 
annealed approximation and the GFs approach are introduced. Presentation of the results is preceded by 
sectionj^where a simpler system-a RC adsorbing onto a homogeneous surface-is considered to illustrate 
the main features of the annealed approximation. This system also serves us as a reference. Section]^ 
is the main part of the paper devoted to the exposition of the results in the case of symmetric compo¬ 
sition and symmetric “interaction map” (“interaction matrix”). We present temperature dependences of 
the total and partial adsorbed fractions, analyze the effect of a variation in the composition and in the 
character of correlations in the RC and the RS, discuss the effect of asymmetry in the composition and/or 
interaction matrix on the run of these dependences. Finally, we summarize in section]^ 


2. Model and method 

2.1. Definition of the model 

Consider a RC chain composed of A and B monomer units interacting with a RS that carries a and 
b sites, figure (note that the two-species heteropolymer chain is similar to the so-called hydrophobic- 
polar (HP) copolymer, a model widely used in theoretical studies of protein folding (43][44l). As in (31), 
we model polymer conformations by 2D-PDWs on the square lattice, therefore, the adsorbing surface is 
simply a line. The monomer sequence of the RC is given hy %- {Xi'X 2 .---.Xn}. where Xi = A or B while 
the surface pattern is given by cr = {...,(7 i,( 72,(T3, ...}, where Ox - aor b. 



Figure 1. 2D partially directed walk model of random copolymer consisting of A (black) and B (white) 
monomer units near the linear (“planar”) random surface composed of a (gray) and b (white) sites. 

Monomer units and surface sites are distributed randomly and modelled as the first order Markov 
chain. For the RC, the latter is determined by the probabilities to find A and B units in the sequence: 
P(Xi - A) - /a and P{Xi - B] - fs - 1- /a, and by the probabilities that the monomer of the type i 
is followed by the monomer of the type j, PiXm - j\Xm-i - i) - Pip The correlation parameter Cp = 
1 “ Pab - Pba - Paa + Pbb - 1 determines the character of correlations in the sequence: Cp > 0 means 
that there is a tendency in the sequence for grouping similar monomers Into clusters, Cp < 0 favors the 
alternating sequence of A’s and B’s, Cp = 0 corresponds to uncorrelated (BernouUian) sequences. The 
probability of occurrence of a particular realization of the sequence x is given by the product P(j) = 

fxiPXlX2PX2X3 '■■■'PXn-lXn- 

Similarly, the sequence of surface sites is determined by the probabilities to find a and b site on the 
surface P{ax - a) - and P((Jx - b] - gh - 1 - gn and the probabilities that the site of the type i is 
followed by the site of the type j, P{(Tx - j\(Xx-i - i) - Sij. The correlation parameter for the surface Cs = 
1 “ ^ab “ Sha - Saa + Sfofo- 1. The probability of a particular realization of the sequence a - {cti,CT 2 , . ■.,cr(} 
is given by the product P{a) = ga^ Sa^a 2 ' ■ • ■' ^oe-iae ■ 

The Hamiltonian (the energy) of the system can be written as follows: 

N 

( 2 . 1 ) 

i=l 
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and the partition function of the systent for particular realizations of the RC and the RS is given by 


■Zn(;6ll,o') = ^exp 
r/ 


N 


i=\ 


( 2 . 2 ) 


where r, = {xi, y/}, i - denotes the chain conformation (the trajectory of the chain), Sij is the 

Kronecker delta, Cij is the energy of the monomer-surface contact (/ e {A, B} and j e {a, b]), and p = 
IIk^T is the inverse temperature. 


2.2. Annealed approximation and generating functions approach 

We will solve the problem in the annealed approximation, where the quenched free energy ppq - 
-((lnZn{P\X:0'))x)a obtained by averaging the logarithm of the partition function over all possible re¬ 
alizations of the RC and the RS is approximated by the annealed free energy = -ln((Z„(;S|x,CT)>;[;>ff, 
where the partition function is averaged prior to taking the logarithm. Here, the angular brackets (...) 
denote averaging over sequence or surface randomness. 

To calculate the annealed partition function of the system considered, {{Zn{p\x<(x))x)a, we use the 
generating functions (GFs) approach (or the grand canonical approach). In the case of adsorption onto 
heterogeneous surface, this approach consists in calculating the GF 

oo n 

H(z,t)=^ ^ ((Z„,„(^ll,CT)>;f>^z”t'”, (2.3) 

^=1 m=l 


where Zn,mip\X’ o') is the constrained partition function of a chain with n monomer units and the length 
of the chain projection onto adsorbing substrate equal to m. The GF variables z and t conjugate to the 
chain length and the chain projection, respectively. The smallest singularity zdt] of Eiz, t] calculated at 
t - \ gives an asymptotic expression for the partition function in the large n limit: ((Z„(/l|x,(7)>^>cr == 
Zp ”(f = 1). Then, the monomer chemical potential (the free energy per monomer unit) p = ln[zc(t = 1)] 
As it was shown in im, the smallest singularity of E{z, t) is associated with the smallest root of the 
equation 

det[E-E;L(z,f)E;s(zf)] = 0. (2.4) 

In equation |2.4| , the functions Hs(zt), !El(z, t) are the GFs of adsorbed segments (usually called “trains”) 
and loops, respectively, in the matrix form: 

OO 

H;s(zr)= ^ns(«)(zr)"R”, 

n=l 
oo n 

“lIz, r)=^ (2.5) 

n=2 m=2 


where Qdn, m) is the number of loops of contour length n and projection length m and klsin, m) is the 
number of trains of length n. The matrices P and S are the transition probability matrices for RC and RS, 
respectively: 


P 

S 


[ Paa Pab ] ^ [ Paa 1 - Paa ] 
i Pba Pbb ) i 1 - Pbb Pbb ) ’ 

I ^aa ^ab 1 ^ ~ 1 

V ^ba ^bb ) V 1 ~ ^bb ^bb ) 


( 2 . 6 ) 


P i8i S is the Kronecker product of the these matrices. The matrix R is defined as 


R= (Pi8>S) W 


(2.7) 


with the diagonal “interaction matrix” 


W = diagiWAa, WAb, WBa, Wsb) 


(2.8) 
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containing all statistical weights of different monomer-surface contacts, Wij = where i e {A, B] 

and j e {a, b]. 

The matrix GFs Es(zt) and El(z, t] can be easily calculated by using the eigenvalues and eigenvectors 
of the matrices P, S, and R and the expressions for scalar GFs of trains calculated straightforwardly as 
follows: 

oo ^2 

“sUl = Y. = z^ + z^ + ---= - -, (2.9) 

«=1 1 - z 

and loops (calculated by using a loop decomposition described in (31)) 


Hl(z, f) 


OO n 


Y Y ^L{n,m)z'^ 

n=2 m=2 


tm-2 


I- zt- z^ - z^t- \/ [I - zt - z^ - z^ -Az'^t^ 
2z^ 


( 2 . 10 ) 


see also (^, equations (23)-(25) for details. 

Note that the determinant equation l |2.4| is a generalization of the analogous scalar equation for a 
homopolymer adsorption onto a homogeneous surface l45l : 1 - E.i^{z]Es{wz) - 0, where w - is the 
statistical weight of a monomer-surface contact, Hl(z) = Hl(z, 1). 

The smallest singularity Zc must then be compared with the smallest singularity zy of the GF for the 
free (desorbed) chain in a bulk H!v(z) = which does not depend on monomer sequence. 

Here, Ov(n) is the number of conformation that a chain of n monomer units can acquire in the bulk. For 
the 2D-PDW polymer model zy = 

In the adsorption transition point Zc = zy, hence, the equation 


det [E - ^LlZy, OHslZy f)] (=1 = 0 


( 2 . 11 ) 


determines the position of the transition point. 


2.3. Calculation of observables 


Various observables can be found via differentiation of the smallest singularity zdt) of the GF. For 
example, logarithmic derivative of Zc with respect to the statistical weight of different monomer-surface 
contacts, ivij, gives the average fraction of these contacts occurring in the adsorbed RC chain: 

dlnzr Wij dzr 

Bi i =-^ =- - -^ . ( 2 . 12 ) 

dlnwij Zc dwij 


The total adsorbed fraction is given by the sum of four contributions: 9 - Bau + BAb + Bbu + Bsb- 

The partial derivative dzcldwjj in equation l|2.12| can be calculated by differentiating the equation 
(2.4|. Namely, if we denote in the left-hand side of equation (2.4| as D := det[E - HlIz, f)Hs(zr)], then 


dzc dDIdwij 

dwij dDIdz ^=zc,t=i 

The derivatives of D can be found with the aid of Jacobi’s formula 


ddetX(T) 

dr 


= Tr Adj(X)- 


dX) 

dr J ’ 


(2.13) 


(2.14) 


where “Tr” stands for the trace of a matrix and Adj(X) denotes the adjugate matrix for X. 

By analogy, one can calculate other observables by choosing a proper differentiation variable. Taking 
a derivative of zdt) with respect to t gives access to the ratio of the average projection of the RC chain 
onto the surface to the RC contour length: 


(m) 

n 


dlnzr 


din t 


t=i 


1 

Zc 


dZc 

dt 


f=i 


(2.15) 


Since in the annealed system, RC monomer units and RS sites are in thermal motion, it is especially 
interesting to find equilibrium moments of their distributions. To calculate the fraction of A monomers 
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in the RC chain, Vyi, let us introduce an auxiliary term into the Hamiltonian (13: we set H ^ H+ AH, 
where 

N 

AH=-(77//3)^Xi. (2.16) 

Later 77 will be set to zero. With this auxiliary term, the matrix P modifies as follows: 


p ^[ Paas^ Pab ] 
i Pbas^ Pbb }' 


(2.17) 


The matrix R is again R = (P isi S) ■ W. Then, the average fraction of A monomers in the RC chain is given 
by 


= - 


dlnzc 


dr] 


r]=0 


1 

Zc 



(2.18) 


The fraction of B units follows automatically: vb-I-va- 

Similarly, one can calculate the average fraction of AA, AB, BA, or BB dyads in the sequence. For 
example, vaa is calculated by introducing the following auxiliary term: 


N 


i=2 


This modifies the matrix P in the following way: 

Paa^"^ Pab ] 
i Pba Pbb j’ 


(2.19) 


( 2 . 20 ) 


i.e., the multiplier appears at the corresponding element (paa) of the probability matrix. The cluster 
parameter Ap, which is the a posteriori analogue of the a priori correlation parameter Cp, is calculated as 
Ap = 1-vabIva-vbaIvb = vaaIva-vbbIvb-1- 

In a similar manner, we can use the same idea to obtain the composition of RS. Thus, to calculate the 
fraction of a sites on the surface, the matrix S should be modified as foUows: 


^ f Sab ^ 

V ^ba^^ ^bb ) 


( 2 . 21 ) 


However, there is an important difference between RC and RS: while in the former case all monomer 
units can be involved in the interaction with the surface, in the latter case only a part of the RS may be 
involved in the interaction with the polymer chain, the size of this “contact zone” is equal to the projection 
of the RC on the substrate m^n. Mathematically, this difference is expressed in the limits of summation 
in the GF of equation In this sense, the quantity “average fraction of a sites on the surface” is not well 

defined for the whole surface because if one takes a finite but large surface which cannot be completely 
covered by the RC there will be two parts of the surface, i.e., a part involved and a part not involved into 
interaction with the RC. A better quantity is the “local” fraction of a-sites, i.e., the ratio of the number of 
a sites in the “contact zone”, nia to the (instantaneous) size of this zone m. 

An expression analogous to equation ^2.18^ does not give the sought fraction of a sites on the surface. 
It gives the ratio of the average number of a sites on the surface “occupied” by the RC, ma, to the contour 
length i.e. nia/ n. Therefore, a correct estimate of the “local” fraction of a-sites, Va, will be given by 


dlnzc ((^l ^ 
dp r,=o\ n j 


Zc drj n ) 


( 2 . 22 ) 


and similarly for the number of dyads. The cluster parameter for the RS is calculated as As = 1 - v^fe/va - 
Vba/Vb- 
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3. Reference system: RC adsorbed on homogeneous substrate 


In the Introduction we briefly discussed the merits and demerits and the fields of application of the 
annealed approximation in the study of (various) disordered systems. It will be instructive to illustrate 
some peculiarities of the annealed approximation in the case of a simpler system, where a RC is adsorbed 
onto a homogeneous surface. This problem was studied earlier: while in (29], the Bernoullian RC was con¬ 
sidered, in (41] a general program for correlated Markovian RC in the framework of the GF approach on 
the lattice was developed. The results in (21 were obtained for a simple random walk model of polymer 
in three dimensions, although the approach is quite general and may be straightforwardly extended to 
other types of lattice models of polymers or to other geometries of adsorbing substrates. In the present 
work, we apply this program for the 2D-PDW model used in the present study. This will also serve as a 
“reference system” for comparison with our “original” system. 

As it was shown m, for a RC adsorbed onto a strictly homogeneous surface, the smallest singu¬ 
larity of the GF E{z) - T.^^i(Zn(l5\x))x^" annealed approximation is found from the equation 

similar to equation (2.4| : det[E- E;l(z)!Es(z)] = 0, with matrix GFs Hl(z) = and Hs(z) = 

^sin)z" R”, R = P-W, where W = diag{WA, wb) - diag(e“^'^'^, is the diagonal matrix of statis¬ 
tical weights of A and B contacts with the surface. Calculation of observables is similar to that described 
is section |231 For more details see (411 . 

These formulas can also be directly obtained from equations (2.3| - (2.8| by assuming that the RS con¬ 
tains only a sites and, correspondingly, by replacing the matrix S in equations by 1 X 1 unity 

matrix, S = 1, redefining eao = £a and cga = £b, and setting r = 1. Then, the interaction matrices W, 
equation (2.8] , and R, equation wfll have the dimensionality 2x2. Alternatively, one can consider a 
and b sites as identical, i.e., set CAa - ^Ab - £a and = ^Bb - ^b and then directly use the formahsm 
introduced in section!^ 

We consider the case of adsorbing A and neutral B contacts: = -1, = 0. Without loss of generality, 

let us assume that A and B monomers have equal probabilities to be found in the monomer sequence, 
fg- 0.5 but have various character of correlations in the RC chain (i.e., various Cp). For the sake of 
comparison, we also consider the homopolymer (HP) consisting of adsorbing A monomer units only. 



(d) 


(e) 


Figure 2. (Color online) Fraction of A- (good) (a) and B- (bad) (b) contacts, total adsorbed fraction (c), 
average fraction of A-monomers, and (d) cluster parameter for RC (e) in the reference system as functions 
of inverse temperature /3 calculated for /a = 0.5 and various Cp as indicated. 
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Figure shows temperature dependences of contact fractions (a, b), total adsorbed amount (c), frac¬ 
tion of A monomers (d), and cluster parameter Ap (e). At Cp > -1, with an increasing p, the number 
of adsorbed A monomers as well as the overall adsorbed fraction monotonously increase while the frac¬ 
tion of adsorbed B units behaves non-monotonously and vanishes at very high p. This is accompanied by 
interconversion of repelling monomer units B into adsorbing units A (however, at ;6 ^ /Itr. where Ptr is 
the transition point, and Ap = Cp) which explains the observed depencences of 9a and 6b- Only 

in the case Cp = -1, corresponding to the regularly alternating AJ3-copolymer its regular (i.e., quenched) 
sequence does not change and both 6a and 9b grow with increasing p due to a cooperative effect. There¬ 
fore, the annealed approximation corresponds to the physical situation that essentially differs from the 
quenched case. 


4. Results and discussion 

4.1. Choosing the system parameters 

The system we consider depends on a large number of variables (to be precise, there are 9 parameters: 
/a. ga, Cp, Cs, CAa, CAby CBay Csb, and p). In IM1 we have suggested a reasonable choice of parameters by 
keeping the inverse temperature p - 1/k^T as a separate control variable and by fixing attractive and 
repulsive energies, thus obtaining physically relevant temperature dependences. In (31), three different 
“interaction schemes” for the monomer-site interaction energies were considered, the choice was also 
motivated by the works of other authors. 

Since the present work is devoted to the study of a particular case of the RC and RS symmetric with 
respect both to composition, fA-ga- 0.5, and monomer-site interactions, CAa - CBb and CAb = eba, this 
restricts even more the set of variable parameters; in fact, there remain only three ones: Cp, Cs, and p. 
Therefore, the main part of the present work will be devoted to the study of the symmetric (but highly 
non-trivial) case; the effect of the asymmetry on the RC and the RS adsorption behaviour will be briefly 
discussed in the end of this section. We will study the system with the following set of parameters: Aa and 
Bb contacts are favorable, CAa - CBb - - l,Ab and Ba contacts are neutral, = cba = 0- 


4.2. Phase diagram and symmetry properties 


In (^, phase diagram, i.e., the dependence of 
the inverse adsorption transition temperature /Itr 
on the correlation parameters Cp and Cs, was cal¬ 
culated. Figure shows this diagram as a density 
plot in (Cp, Cg) coordinates (in such a form it was 
not presented in (STj). As it follows from the di¬ 
agram, the smallest values of Pti (dark color, the 
bottom of the color scale in figure]^ are observed 
in the vicinity of (Cp, Cg) = (1,1) and (-1,-1) 
whereas the largest values of jStr (light color, the 
top of the color scale in figure are observed 
in the vicinity of (Cp, Cg) ~ (-1,1) and (1, -1). We 
can also observe that this density plot is symmet¬ 
ric with respect to the origin (Cp, Cg) = (0, 0). This 
means that PtriCp = x, Cg = y) = PtiiCp = -x, Cg = 
-y) where x and y may take on any value in the 
interval (-1,1). This symmetry of the diagram is 
the consequence of the interaction and composi¬ 
tion symmetries. (At a glance it may also seem that 
the diagram is also symmetric with respect to the 
diagonals Cg = +Cp but this is not the case). 





Figure 3. (Color online) Inverse transition tempera¬ 
ture in the annealed approximation as function of 
RC and RS correlation parameters for - ^Ba - 
-1, CAb - CBa = 0 in the symmetric case fA = ga = 
0.5. 
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This symmetry manifests itself in temperature dependences of the observables beyond the transition 
point. We vary the values for the correlation parameter for the RC and the RS considering the cases 
Cp = -0.5, 0, O.S, and Cs = -1, -0.5, 0, 0.5. It is clear that various combinations of Cp and Cs are possible 
(3 X 4 = 12 combinations). However, due to the interaction and the composition symmetries, it turns out 
that these 12 dependences for an adsorbed fraction can be presented on the same plot by 5 curves. 


4.3. Adsorbed fraction 

Figure (a)-{c) shows the dependences of the fraction of good (Aa, Bb) and bad (Ab, Ba) contacts 
and of the total adsorbed fraction, respectively, on the inverse temperature /3. Due to the symmetry of 
the system, we have Oau - dsb and dAb - Gsa- With increasing (5, the fraction of good contacts grows, 
the maximum “saturation” value for both Aa and Bb contacts equals O.S, which is the maximum that is 
available a priori (according to the values of /a and ga), i.e., that A (5) monomer units and a (b) surface 
sites may form. For the total amount of good contacts 0good = ^Aa + dsb^ the upper boundary Is equal to 
unity. In other words, this indicates that A^ B and a^b transformations do not occur, as opposed to 
the reference system where the fraction of good A-contacts grows as p increases, figure[^(a). The fraction 
of bad contacts, figure (b), behaves non-monotonously: just above the adsorption transition it grows 
with an increasing p but then decays to zero. Since the fraction of bad contacts is much lower than the 
fraction of good contacts, the overall adsorbed fraction increases monotonously with p, figure|^{c), and 
the energy per monomer unit, figure [^(d), monotonously decreases from 0 to -1. All the curves start In 




(a) (b) 




(0 


(d) 


Figure 4. (Color online) Fraction of good (a) and bad (b) contacts, total adsorbed fraction (c), and energy 
(d) as functions of inverse temperature p calculated for /a = ga = 0.5, and various values of cluster 
parameters-curve 1: Cp = -0.5, Cs = -1; curve 2: Cp = +0.5, Cs = +0.5; curve 3: Cp = 0 or/and Cs = 0; curve 
4: Cp = +0.5, Cs = +0.5; curve 5: Cp = 0.5, Cs = -1. 

the transition points ranged in accordance with the phase diagram, figurej^ 

4.4. Moments of monomers' and sites' distributions 

In order to understand the behavior of the observables (which is, as we see, much more reasonable 
and closer to what one can expect In the quenched system, as compared to the reference system) bet- 
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ter, we study the temperature dependences of a posteriori moments of distributions of monomers and 
sites. The first remarkable result is as follows: in the considered symmetric case, the fractions of A (B) 
monomers and a (b) sites do not depend on the temperature and are always equal to O.S! This takes place 
at any combinations of correlation parameters Cp and Cg (therefore, we do not show these dependences 
in the figure due to their trivial form). 

As regards the second moments of the a posteriori distributions, these are changing with 15. This 
means that there are transformations in the RC and the RS sequences but since the RC and the RS compo¬ 
sitions are invariant, the transformations occur according to a specific law: A and B units (a and b sites) 
move in the RC (RS) in order to tune their sequences with respect to each other in the best way to reach 
the lowest interaction energy. Alternatively, these rearrangements can be considered as coupled chem¬ 
ical reactions: for example, a transformation (Xi = A) ixi - B] should occur simultaneously with the 
reaction {%] -B)^ ixj - A), i and j + i denote positions of the monomer unit in the RC sequence. At the 
same time, another important rule sill holds: in the transition point, the first and the second moments of 
monomers and sites distributions are equal to the corresponding a priori probabilities. 

Now, let us consider in detail some particular cases of RC adsorption onto RS. 


4.4.1. Quasi-blocky RC 

As it follows from figures[^and[^ a quasi-blocky RC has a better capability of adsorbing onto a quasi- 
blocky RS rather than onto a quasi-alternating RS. As p increases, the tuning of the character of corre¬ 
lations in both RC and RS occurs, figure]^ (a), (d): When Cg - 0.5 (= Cp), both Ap and As slightly increase 
whereas at Cg - -0.5 (= -Cp), both Ap and As move towards each other, i.e., the positive Ap decreases 
and the negative As increases. At Cs = 0, an interesting behavior is observed: the RC cluster parameter 
does not change but the RS cluster parameter increases. At Cs = -1, the surface regularly alternates (i.e., 
quenches) which cannot change “by definition”, and the RC should adjust itself; hence. As decreases and 
becomes negative. 



(a) (b) (c) 





p P ft 

(d) (e) (f) 


Figure 5. (Color online) Cluster parameters for RC (a)-(c) and RS (d)-(f) as functions of inverse tempera¬ 
ture ft calculated for fA = ga = 0.5, Cp = 0.5 (a, d), -0.5 (b, e), 0 (c, f) and various Cg as indicated. 


4.4.2. Quasi-alternating RC 

Here, the same tendency to tuning as in the previous case is observed, figure (c), (f). Hence, all 
the corresponding arguments may be straightforwardly reproduced with “the change of the sign”: at 
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Cs = -0.5 (= Cp) both Ap and Ag slightly decrease, at Cg = 0.5 (= -Cp) Ap and Ag move towards each other, 
at Cs = 0 Ap does not change and As decreases, at Cs = -1, Ap decreases towards -1. 

4.4.3. Bernoullian RC 

With respect to the transition point and the adsorbed amount, figures]^ andthe behavior of the 
annealed Bernoullian copolymer does not depend on the correlations in the annealed RS (the latter can 
be even quenched regularly alternating with Cs = -1). For the second moments of monomers and sites 
distributions in the RC and RS, respectively, the following regularity is observed: the a posteriori cluster 
parameters As in the RS — the “partner” of the Bernoullian RC — does not change with p and remains 
equal to the a priori correlation parameter Cs, figure[^(b), (e). The "double Bernouhian” case Cp = Cs = 0 
is remarkable: the second moment does not change both in RC and RS. 

4.4.4. Summary for the symmetric case 

The analysis of the first and the second moments of monomers and site distributions in RC and RS, 
respectively, in the symmetric case shows that upon a decrease in temperature (increase in p), the in¬ 
crease in the number of good contacts and the corresponding decrease in the number of bad contacts is 
implemented via rearrangement of ^ and B monomer units (a and b sites) in the RC chain (on the RS). We 
have also seen that the first moments (and, in some cases, the second moments) keep their expected val¬ 
ues and, in this sense, the annealed approximation turns out to be more accurate than in the case of the 
(more simple) reference system. This also means that an improvement of the annealed approximation 
with the aid of the first-order Morita approximation (and the second order-in the case of the Bernoullian 
RC and RS) will produce no effect because the corresponding Morita constraints are already satisfied in 
the annealed system and the same results wih be obtained. 

4.5. The effect of asymmetry 

As we have seen, in the symmetric case, the adsorption of annealed RC onto annealed RS exhibits the 
most interesting features, the most essential among them being the correct (with respect to the a priori 
probabhities) values of the first moments of monomers and the distributions of sites meaning that the 
overall composition of the RC and RS does not change during the interaction. This is the consequence of 
the system symmetry. When both the RC and the RS are, in addition, Bernoullian, the second moments 
keep their correct values too. 

Let us now analyze how a violation of the system’s symmetry affects the RC adsorption onto the RS. 
In principle, the symmetry can be destroyed in two different ways: one can either (1) choose non-equal 
probabhities for A and B units (a and b sites) to appear in the RC sequence (on the RS), i.e., set /a, ga + 0.5 
or (2) change the value of one of the interaction parameters e/y. To see how different ways of introducing 
the asymmetry influence the adsorption, let us take the (simplest) set of parameters with /a - g^- 0.5, 
Cp = Cs = 0 (both RC and RS are Bernoullian), CAa - ^Bb - “1. and CAb - ^Ba = 0 as the reference. We whl 
change one of the composition or the interaction parameters, keeping other parameters unchanged, as 
in the reference system. In particular, we (1) slightly change the probabilities of the appearance of A and 
B monomers in the RC, /a - 0.6, keeping for the RS ga = 0.5 or (2) increase the attraction of A monomer 
to a site by setting CAa = -1.2 or (3) make one of the bad contacts repulsive instead of neutral by setting 
CAb = -t 1, or (4) turn the good Bb contacts into bad by setting csb- ^Ab - ^Ba - 0. 

These four cases are compared in figure]^ The results of the comparison are as follows: (1) an increase 
of the probability of A monomers present in the sequence, changes neither the position of the adsorption 
transition point nor the shape of the temperature dependence of the adsorbed fraction (and of the GF 
smallest singularity). At the same time, the RC and the RS compositions behave differently with the change 
of temperature: as the inverse temperature p increases, composition of the RC remains the same whereas 
the fraction of a-sites increases and tends to 0.6. (2) An increase in the Aa affinity results in the shift of 
the transition point toward smaller p, while with an increase of p, the (expected) growth is observed in 
the amount of A monomers and a sites providing the most favorable Aa and Bb contacts. (3) An increase 
in the Ab interaction energy, shifts the transition point to larger values of p-, fractions of A monomers 
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Figure 6. (Color online) Overall adsorbed fraction (a), fraction of A-monomers (b) and a-sites (b) as 
functions of inverse temperature p calculated for symmetric reference system (red solid curves) with 
fA = Sa- b-5> cp = Cs = 0, - “b and = 0. Other curves are for the systems differing 

from the reference one by one of the parameters, indicated at the curve. 




(a) (b) 


Figure 7. (Color online) Overall adsorbed fraction (a) and fractions of A-monomers (solid curves) and a- 
sites (dashed curves) (b) as functions of inverse temperature p calculated for Cp = Cg = 0, - ^Bb - “b 

£Ab = eBa = 0, and equal fractions of A-monomers and a-sites (.fj{ = ga) as indicated. 


and B sites slightly decrease just after the adsorption transition; as p increases, they both tend to 0.5, in 
order to have a maximum of possible favorable Aa and Bb contacts. (4) The transformation of attractive 
Bb contacts intro neutral ones shifts the transition point towards larger p and leads to massive B ^ A 
and a transformations with an increasing p. 

The observed behavior is in agreement with the behavior of the reference system (RC adsorbing onto 
a homogeneous surface, see section [^, where the equilibrium is shifted towards transformation of non¬ 
adsorbing monomer units into adsorbing ones. In our case, we encounter more comphcated coherent 
monomers and transformations of sites. As soon as some preferences appear in the transformation re¬ 
action constants (governed by /a and ga) or in the interaction map, the equhibrium shifts toward these 
preferences. 

There is another interesting way to partly violate the symmetry of the system. By comparing the 
fraction of A-monomers with that of a-sites and the fraction of B-monomers with that of b-sites we can 
obtain the upper boundary for the maximum possible fraction of good contacts. Since good contacts are 
the Aa and the Bb ones, then 0good = Sao + dBb- Obviously, dAa ^ nrin(/A, ga) ydBb"^ min(/B> hence, 

= min(/A,ga) +min(/B,gb). In the symmetric case, = ga = gs = 0.5, therefore, = 1. If 

we choose the RC and the RS compositions, so that /a - ga and fs - gb, the upper boundary for the total 
fraction of good contacts whl still be equal to one, as it is in the symmetric case. 

Temperature dependences of the overall adsorbed fraction, RC and RS compositions and cluster pa¬ 
rameters for /a = ga s 0.5 and Cp = Cg = 0 are presented in hgure[^ We see that in spite of the invariance 
of 0™^, a simultaneous increase in the fraction of A-monomers and a-sites favors adsorption: the tran¬ 
sition point shifts to the lower p values, for a given temperature, the total adsorbed fraction is larger in 
the case of larger /a - ga- Moreover, with an increasing p, the fractions of the major A and a components 
grow and the fractions of the minor B and b components correspondingly decrease. It is also clear that if 
we take /a = ga < 0.5, the picture will remain qualitatively and quantitatively the same with respect to 
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major and minor components (here, B and b are major components, A and a are minor components). The 
oniy exception is the symmetric case fj^- g^- 0.5 , where there is a strict A/B and a/b balance. That is, 
one can say that the equilibrium gets shifted towards the major component, and the situation becomes 
similar in some sense to that in the reference system of sectionj^ 


5. Conclusion 

We have considered a two-dimensional partially directed walk (2D-PDW) model of a two-letter {AB) 
random copolymer (RC) adsorption onto a two-letter (ab) random surface (RS). This model was introduced 
in our previous work IMl where it was treated by using the combination of the annealed approximation 
(to perform double averaging over the sequence and surface disorder) with the generating functions 
(GFs) approach (to sum over polymer conformations). In contrast to (31), in the present work we have 
gone beyond the transition point and studied the temperature dependences of various observables for the 
annealed symmetric system. This choice was motivated for several reasons: The annealed approximation 
provides a zero-order rough approximation to a reahstic quenched system; on the other hand, it can serve 
as a prototype of real physical systems like two-state polymers (^. The system was chosen to be “two¬ 
fold” symmetric: with respect to the composition of the RC and RS and with respect to the interaction 
map (good attractive Aa and Bb contacts had the same energy, while both bad Ab and Ba contacts were 
equally neutral). In (31) it was shown that in this case the system has the most interesting phase diagram. 
Finally, the symmetry of the system makes highly unpredictable interconversion of monomers and sites 
(A^B and a^b)m the annealed system. Therefore, special attention was paid to a posteriori moments 
of distributions of monomers and sites. 

We have shown that in the considered symmetric case, the expected a posteriori compositions of the 
RC and the RS correspond to the a priori probabilities to meet A monomers and a sites in the RC monomer 
sequence / RS site sequence and do not change with the temperature. At the same time, the a posteriori 
cluster parameter (related to the second moments of distributions of monomers and sites) in the RC and 
RS changes with the temperature, indicating that monomers and sites rearrange in the RC and the RS 
to provide a better matching between them and, hence, a stronger adsorption. A special case is the one 
where both the RC and the RS are Bernoullian: in this situation, both first and second moments keep their 
correct values at any temperature. 

We have also studied the effect of the system symmetry violation on the adsorption behavior. There 
are various ways of doing this and all of them shift the equilibrium towards the major component and/or 
more favorable contacts. 
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AACopGi^m ci/iMeTpi/iMHoro Bi/inaAKOBoro Kono/iiiviepa Ha 
ci/iMeTpi/iHHy Bi/inaAKOBy noBepxHio: BiAna/ieHi/iii Bi/inaAOK 

A.A. no;ioLibKHbP2 

^ iHCTi/iTyr MaKpoMO/ieKy/inpHi/ix cnonyK, PociPicbKa aKafleMin nayK, 199004 CaHKT-nerepOypr, PociiiCbKa 
tPeflepapin 

^ CaHKT-fleTepOypsbKi/ifi HapioHa/ibHi/iit flocniAHupbKHM yHiBepci/iTer iHc|)opMai4iMHi/ix rexHO/iorifi, MexaniKi/i ra 
oniMKi/i (yniBepcuTeT ITMO), 197101 CaHKT-neTep6ypr, PociilcbKa (UeAepaLtin 

AflCopOpin CMMerpi/iHHOro {AB) Bi/inaflKOBOro Kononiwepa (BK) na ci/iMerpnnHy {ab) Bi/inaflKOBy HeoflHopiflHy 
noBepxHKD (Bfl) Bi/iBnacTbcn y Had/inxtenHi Biflnany is Bi/iKopucraHHnMi/i fl/in noniwepa flBOBi/iMipHOi MOfleni nac- 
TKOBO HanpnM/ieHnx 6;iyKaHb. flOKasaHO, tpo y CMMerpi/iHHOMy Bi/inaflKy OHiKyBani a posteriori KOHiteHTpapii 
BK i Bfl MaiOTb npaBi/iabni SHaneHnn (nxi BiflnoBiflaiOTb Tx a priori MMOBipnocTriM) i He SMiHioiOTbcn, b saaextHO- 
CTi Bifl TeMneparypi/i, b tom nac, nx flpyri MOMeHTi/i posnofliaiB MOHOMepiB i BysaiB b BK i Bfl SMiHioiOTbCfl. Lfe 
noxasye, 040 MOHOMepw i Bysai/i BsacMHO He nepeTBoproHDibcn, a ni/itue neperpynoByHDTbcn, 0406 saOesneni/iTM 
xpaute AonacoByBaHHfl Mi>x hmmi/i i, nx pesynbiar, ci/inbHiujy aflCopOpiK) BK na Bfl. flpore Oy^b-nxe nopytueHHn 
CMMeTpii CMCTeMM scyBac piBHOBary y HanpnMxy ochobhoi xoMnoneHTM i/a6o /40 6mbuj cnpHnTtii/iBi/ix xonraxiiB 
ra npi/iBO/41/iTb flo BsaeMonepeiBopeHHn woHOMepiB i BysniB. 

KnioMOBi cjioBa: Bi/inaflKOBi/iM KononiMop, Bt^naflKOBa noBopxHB, noniMepna aflcopBpin, Ha6ni/i>KeHHn 
Bifinany, renepyHMi cpyHKpii 
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